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Abstract. We study the Drude weight D(T) at finite temperatures T of 
an integrable bosonic model where the particles interact via nearest-neighbour 
coupling on a chain. At low temperatures, D{T) is shown to be universal in the 
sense that this region is equivalently described by a Gaussian model. This low- 
temperature limit is also relevant for the integrable one-dimensional Bose gas. We 
then use the thermodynamic Bethe ansatz to confirm the low-temperature result, 
to obtain the high temperature limit of D(T) and to calculate D(T) numerically. 



PACS numbers: 05.30.-d, 05.30.Jp 

1. Introduction 

The one-dimensional Bose gas with J-interaction was solved exactly by Lieb and 
Liniger [1] ; thermodynamic properties were calculated by Yang and Yang [2] , using a 
method which is nowadays known as thermodynamic Bethe ansatz (TBA) . This Bose 
gas model can be viewed as the continuum limit of a one-dimensional lattice model, 
where the particles interact via a short-range interaction. The advantage of a lattice 
description is that integrals over distances (momenta) have a natural lower (higher) 
cutoff, namely the lattice constant (bandwidth). This makes it possible to calculate 
high-energy properties, in our case, the high-temperature limit of D{T). In the other 
extreme, at low energies, a continuum description of the lattice model is justified, 
which opens the possibility of covering the low-energy regime of the continuum model 
as well. 

Here, we want to benefit from the existence of an integrable lattice regularization 
of the Lieb-Liniger (LL) Bose gas, stemming from a g-deformation of the underlying 
bosonic commutators. This (/-Bose model was shown to be integrable in [3, 4, 5], where 
the close relation to the Lieb-Liniger model in the continuum limit was also pointed 
out. Very recently [6], the lattice model has been rederived in the framework of the 
algebraic Bethe ansatz. In that approach, the particle current as a conserved quantity 
appeared naturally. 

In this work, we study the dynamical conductivity cr(w,T) of the q-Bose model, 
depending on the frequency uj and the temperature T. It is defined by 

J{uj,T)^a{uj,T)E{io), (1) 

where is the particle current induced by the field E. The particle current is 
conserved [6] , leading to 

Re cr(w,T) = 7r(5(a;)D(r), (2) 
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that is, the conductivity is infinite at zero frequency. Our aim here is to calculate the 
temperature dependence of the Drude weight D{T). 

This paper is organized as follows: In section 2, wc shortly review the model and 
provide some background on the calculation of a{uj,T). The third section is devoted 
to the calculation of the leading term in D{T ~ 0) from the effective low-energy 
Gaussian field theory. Section 4 contains the calculation of D{T) from the TBA. 
The T = 0-result of the preceding section is recovered here, as well as an analytical 
formula in the high-temperature region £)(T 3> 1). Numerical results for a wider range 
of temperatures interpolate between these two extremes. We end with a conclusion 
and outlook in section 5. 

2. Conserved current in the g-Bose model 

In this section, wc shortly review the g-Bose model following [6]. Let b, b''^,n be the 
generators of the g-oscillator algebra: 

bb+ = 1 - r?2"+2 ^ btb = 1 - , 

bn=(n-M)b, b^n = (n - l)bt , 
with < g < 1. In the following, it is often convenient to use the parameterization 

q = T] e ]0, oo[ . 
A g-oscillator is represented in the Fock space: 

n\n)=n\n); n = 0, 1, 2, 3, . . . ; b|0) = . (3) 

We deal with a chain of length L where the local g'-oscillator algebra is assigned to 
each site i= 1, . . . ,L. The quantities 

1 ^ 1 ^ 

^ i=i ^ t=i 

commute pairwisely. The eigenvalue of AA is the number of particles N. The choice of 
the Hamiltonian 

n = -l{p+ + r-)-fiM, (4) 

where /x is the chemical potential, leads [3, 4, 5] in the continuum limit to the Lieb- 
Liniger Bose gas with Hamiltonian 

W = ^ y" [d^'^Hx)da^'^{x) + 2c^''f(x)*'f(x)4'(x)4'(x) - dx. (5) 

Here A is the lattice constant of the g-Bose model and the coupling constant c has been 
introduced such that r] = c - A. Furthermore, trivial constants have been absorbed 
into 11. 

The current operator J' is defined by 
dt^e = ji-i - je, 
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where 9tn^ = i \H, n^]. It turns out that 

ji = b^b^^+i - he+ih^e, (6) 
so that J = Yi (^+ ^ ^- ) ' ^'^d thus 

[n,j]^o, (7) 

i.e. ^7 is a conserved quantity. 

Let us finally write the current operator in the same continuum limit that lead 
from (4) to (5): 

j - da;, 

which coincides [7] with the total momentum operator of the model (5). This equality 
of the total current and momentum operator is, however, not true for the lattice model 
[8]. 

Within linear response theory, the current is proportional to the field, where the 
constant of proportionality is the conductivity. It is related to the current-current 
correlation function (Kubo-formula, [9, 10]) as follows: 

/•oo 

a{uj,T)= dte'^' dT(J(-i-ir)J), 
Jo Jo 

where t is the (real) time, (3 := 1/T, and (• • •) denotes the thermal average. Because 

of Eq. (7) we have 

a{u,T)^ lim .. ^ . A J^), 
e^o i(a; - le) 

from which Eq. (2) follows with 

D{T)=p{J^). (8) 
The latter expectation value is obtained straightforwardly from 

{J^)=dl\ntre-^^-^^\^^^. (9) 
As mentioned above, the model defined by the Hamiltonian (4) is Bethe ansatz 
(BA) solvable. Because of Eq. (7), this is also true for 

n:=H- ~XJ, (10) 
where A = A//3. Following the sketch of the BA solution in [6], it is not difficult to 
show that the corresponding eigenvalues E read 

N 

e{k) = — cos k — Xsink — fi, 

where the quantum numbers (BA-roots) fc„ are to be determined from the set of 
coupled algebraic equations 

i,^^ yr sin[(fc,.-fc^)/2 + ir/] 

ll.sin[(fc,-fc„)/2-i^]'^ 

In the thermodynamic limit, ground state and thermodynamic properties of Ti. have 
been determined from these equations [3, 4, 5, 6]. In section 4, the Drude weight 
D{T) will be calculated along similar lines. Let us finally note that the continuum 
limit which lead to the bose gas (5) corresponds to the rational limit k^ — Pn ■ A, 
T] = c ■ A, L = L/A in (11) at fixed particle number N, where L is the length of the 
system. This limit leads to the BA equations found by Lieb and Liniger [1]. 
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(13) 



3. The low- energy effective theory 

In this section, an effective field theory is constructed for the model (10) in order to 
calculate the leading term in a low-temperature expansion of D{T). 

A density-phase representation of the operators b, b'l' has been used in [6] to arrive 
at a Gaussian model which describes the low-energy, i.e. low-temperature properties 
of the model (4). In the field- theoretical calculation we also restrict ourselves to a 
deformation parameter q which is close to 1, i.e. 0. However, the final results will 
be valid for the whole range -q € ]0, oo[. We will come back to this point at the end of 
this section. 

Following [6] I, a continuum limit different from the one leading to (5) is 
performed. Therefore we set a: = £ • A, = p{x) ■ A, where A is the lattice constant. 
Using a density- phase representation of the bosonic operators b^,b^^, one obtains an 
infinite series, of which we only give the leading terms here: 

^ = b(.) = .jm^) (l - ^ + O ((p,A)^)) e-(^) (12) 

^ = bt(.) = e--(^) (l - ^ + O ((p,A)^)^ 

Furthermore, we split p{x) = n + $(a;) where n is the average particle density per unit 
length. iS The density operator p is canonical conjugate to the phase operator 9, such 
that ['^{x), 0{y)\ =\.5{x — y). Since <^{x) has been introduced to encode the low-energy 
density fluctuations, we require {<^{x)) <C n, {dx^{x)) <C 1. To gain insight into the 
physical significance of 0{x), we observe that, following (6), the local current density 
J{x) is written as J{x) = ndx6{x), so that dxO{x) stands for the local momentum 
fluctuations. 

Under rescaling 

0^9' l^/Yc, = ^/kI^' , (14) 

where the charge velocity Vc and the Luttinger parameter Kc are, in this 
approximation, given by (r/ = cA <i; 1) 

(15) 

the Hamiltonian density describing the low-energy fluctuations is written as (we 
discard a trivial additive constant) 

{dx9'f + ($')'] (16) 



Vr 

Hq = — 

2 



H = Ho - py/K^^' - X^=dx0' . 

Let us write n = Vc- Kc in the last term so that the only parameters that enter in the 
Hamiltonian are Vc, Kc and the potentials p, A: 

H = Ho- p\^:<P' - Xvc^/KcdxO'. (17) 
From the Lagrangian Lq corresponding to Hq, the equation of motion 

dl9'/vc - vcdle' = (18) 
X We use a slightly different notation here. 

S In the framework of the BA solution, the same symbol is used for the average particle density per 
lattice site. 
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is obtained, where diQ' jv^ = Eq. (18) can be interpreted as a continuity equation 
for a current v^d^O' and a density dfO' /vc = $' (the existence of a conserved quantity 
with these current and density is due to the invariance of the Lagrangian under 
translations ^ — > 6 + a). This is consistent with our interpretations of 6 and i> above. 

From (17), we can infer the zero-temperature charge susceptibiHty Xc and 
current-current-correlation function {J^). Namely, upon shifting ^ + ^^(^ , 
8x0' dx0' + Xy/K^ we obtain 



H = Ho 



2Vr. 



Thus 



so that 



Xc{T = 0) = ^ , (^')/^ = ^cKc, (19) 



D{T = 0) = VcK^ . (20) 

The Gaussian model (17) is written in the standard form in terms of left- and 
right-moving fields (f>L,R by introducing cp' such that = dx(j>' and (j)' = (j^L + (I^r, 
e' = (t>L-4>R, [11]- Then 



H 



{dx(i)Lf + {dx(i)Rf - fi\/K^dx {(l)L + <Pr) - Xvc\/K'cdx {(I^L - (1>r) ■ (21) 



A few comments arc in order at this point. First of all, Eqs. (19), (20) have been 
derived for 77 ^ 0, resulting in (15). However, as will be shown in the next section, 
Eq. (20) is valid for the whole range rj G ]0,oo[, where Kc have to be calculated 
from the exact solution (the same is true for Xc given in Eq. (19), as shown in [6]). In 
the two special cases r/ — > 0, 77 — > 00, explicit results are available: 

For the latter case, the results for the parameters Vc, as calculated in [6], have 
been inserted. 

Secondly, the results in (19) are not specific for the q-deformcd Bose model 
considered in this paper, but are valid for all models that, in the low-energy region, 
are described by the Gaussian model (17). Let us give two examples. 

The first example is the Licb-Liniger Bose gas (5). As shown by Haldane [12], its 
low-energy properties are equally described by the Gaussian model (17). Especially, 
the relations (15), valid in the small-coupling regime, are the same for the Bose gas (5) 
for weak coupling c. Thus both models arc equivalent in the weak-coupling limit. In 
the strong-coupling regime of the Bose gas, however, Vc and Kc have to be calculated 
from the Bethe ansatz solution [1, 7]. 

As a second example, consider the critical spin-1/2 XXZ-Heisenbcrg chain. 
In this case, the low-energy theory is obtained by performing a Jordan- Wigner 
transformation of the spin operators, resulting in a lattice model of spinless fermions. 
In the continuum limit, the corresponding fcrmionic operators are written in terms of 
bosonic fields (bosonization), resulting in the Hamiltonian density [13, 14] 



lil + {dx<Psf 



TT V TT 



Hxxz = Ho- hJ —dx<ps - M —v,U,, (23) 
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with [(j)s{x),'Rs{y)\ = — y). Here, Vg is the velocity of spin excitations and Kg 
the corresponding Luttinger liquid parameter. In (23), two terms additional to Hq 
have been included: The magnetic field h couples to the spin density (stemming 
from the spin operator Sf in the lattice model), which yields the fermionic density. 
Secondly, A couples to the spin, i.e. fermionic current, originating in the operator 
{S^S^j^-y — S'^ S'^j^-^) / {2i) in the lattice model. By comparing with Eq. (17), one has 
$' = dx(j)' oc dx(i)s and dxO' oc lis. We now set (j)s — 4>l + </>_r, Hs =: dxOs = 
dx {4'L — 4>r)i resulting in essentially the same Hamiltonian density (21) (apart from 
trivial prefactors in the potentials /x, A). Thus the magnetic susceptibility Xs and 
the Drude weight Dg of the spin current are given by Xs{T = 0) oc Kg/vg [13], 
Ds(T = 0) oc VgKg [15, 14], with trivial constants of proportionality. 

We finally remark on next-leading contributions to the leading low- temperature 
result (20). Such contributions are determined by irrelevant operators in the effective 
field theory, which have been neglected in (16). In the case of the g-Bose model, these 
stem from 

i) higher order-terms in the series expansions (12,13), 

ii) corrections to the long-wavelength approximations (12,13) due to the discrete 
nature of the particle density. 

As pointed out in [12, 16], the second kind of corrections can, in lowest order, 
principally give rise to a term oc cos27r((/) + nx/A). Due to the lattice, such a term 
is commensurate at density n = 1. After the scaling (14), we see that this term 
becomes relevant for Kc < 2/tt. The analysis in [6] shows that Kc E ]l/7r,oo[, so that 
a massive phase seems to be possible. Such a scenario has indeed been identified in 
the generic Bose-Hubbard model [17]. In our case, however, this would contradict the 
conservation of the current, Eq. (7): Since [exp [i(/)(a;)] ,J{x)] = S{x — y)nexp[i(j)(x)], 
the above cos-term cannot occur at all.|| This leads us to the conclusion that terms 
due to i) in the above list cancel those due to ii). It would be desirable to show this 
explicitly, a goal which, at this point, is left for future work. 

Further candidates for irrelevant operators have been given in [6]. These have 
even scaling dimensions, so that the Drude weight behaves as 

D{T > 0) = VcKc + O (T^) . (24) 

The amplitude of the next-leading term is unknown so far and given by the amplitudes 
of the leading irrelevant operators. 

It is interesting to note that the cos-terms due to ii) can be "fermionized" and 
are then interpreted as Umklapp-scattering processes. Such processes are present 
in the XXZ-cha.m and accordingly the corresponding operators occur as irrelevant 
contributions in the low-energy effective theory. Additionally, operators due to the 
band curvature at the Fermi points are present. Both types of operators lead to 
temperature-dependent corrections to Dg{T = 0), with non- integer and integer powers, 
respectively [14, 18]. Thus although both models are, in the leading order, equivalently 
described by (21), they differ in the irrelevant operators. 

4. Thermodynamic Bethe Ansatz 

In this section, an exact formula for the temperature dependence of D{T) is presented. 
This formula is evaluated analytically in the low- and high-temperature limits, as well 

II The same is true for higher order terms of the form cos 2iTm(p. 
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4-1. Low- and high-temperature limits 

Employing the TBA originally developed by Yang and Yang to calculate the 
thermodynamics of the Lieb-Liniger Bose gas [2] , the grand-canonical potential of the 
model (4) has been obtained [5] and analyzed [6]. Analogously, the thermodynamic 
potential g{T, jj) corresponding to (10) with A = A//3 reads 

J^A{k) — , (25) 

with a:=l + A and 

J — TT 

sinh 2?7 



Ina(fc) =/3cosA; + /?/i + AsinA;+ / K{k - q)\n A{q)^ (26) 

27r 



cosh 2?7 — cos k 

Then, following Eqs. (8,9), 

D(T)=/3a|/35(r,/z)|,^o_„_,. . (27) 

Note that the derivatives with respect to A are to be taken such that the particle 
density n is held constant. 

Before diving into the analysis of the low- and high-temperature limits, it is useful 
to envision the following relations: 

^ ' 1 -|-exp[-lno(A;)] ^ ' 

^ ^ ' l + exp[-lna(fc)] ^.cosh^'^ ^ ' 

Let us now first evaluate Eq. (27) in the low-temperature limit. In this case, for A = 0, 
lna(A;) = - /3eo(A:) -h O (T) (30) 



K.{k - q)Eo{q)- 

-B 



2-K 

so{±B) = 0, eo{\k\ <B)<0 (31) 
We neglect the next-leading correction 0{T) in Eq. (30) in the ongoing. Then 

K{k — q)d\ In A{q) — 

^-^°sinfc-|- / K{k — q)dxlna{q) — , 
J-B 27r 

where in the last line, we used Eqs. (28,30,31). Obviously, 

Um^ dx In a{k) = dkSo{k). (32) 

Next, we consider the equation for 9^ lna(A;), by making use of Eq. (29): 

dllnaik)^ nik-q)dllnA{q)^ 

' TdkEoiB) [K{k -B) + K{k + B)]+ [ K{k - q)dllna{q)p-. 

J-B 



(33) 
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The last term in the last line is obtained from the first term in (29), using 
Eqs. (30), (31). To obtain the second term, we performed the approximation 

lim K—rrr = Hm ^ „ = ^ \ , \5{k - B) + 6(k + B)] 

4cosh2iil#l /3--4cosh2^ dkSoiB)^^ ' ^ 



where use was made of Eq. (32). One now observes that 

(34) 



/. 



dl\na{q)p-=Tdkeo{B) 

B 



-- + 2p{B) 

TT 



where the function p{k) is given by 

In fact, p{\k\ < B) describes the density of BA-roots in the interval [— B, B]. Eq. (34) 
can be established by writing down the integral series for the left and right hand sides 
of the equation and comparing term by term. Combining (27) with (25) and (33), we 
thus arrive at the result 

D{T>0) = 2p{B)duea{B) (35) 

Now we recall how the charge velocity Vc and the Luttinger parameter Kc are obtained 
within the BA-solution [7, 6]: 

_ dkSo{B) 
27rp(B) 
ire = 4V(B), 

which confirms that Eq. (20) is valid for the entire range rj G ]0, oo[. 

We now consider the high-temperature limit T ^ oo. In this case, Eq. (26) 
simplifies to 

\na{k)=Pp+ f K{k - q) In A{q) — , (36) 

which is solved by 

Ina = -ln(e"'^^^ - 1) , (37) 

without any dependence on the spectral parameter k (note that J^^ K{k — q)dq = 2it). 
In order to derive Eq. (36), it was essential that the fc-dependent driving term in 
Eq. (26) is bounded, a fact due to the presence of the lattice. In the analogous 
equation for the continuous Bose-gas, the driving-term is unbounded [2], which up to 
now makes it impossible to write down a high-temperature expansion for that model. 
Taking into account Eqs. (28,37), we obtain the equation for d\\n.a{k), 



dx\na{k)=smk + el^i' [ K{k - q)dxlnaiq)p-. 

J-7V 277 



(38) 



It is solved explicitly for ^ oo and rj ^ 0: In the rj oo-case, k = 1 
and d\lna{k) = sinfc. In the other extreme, r/ — > 1, K{k) = 27rS{k) and thus 
d\ Ina(fc) = sink/ (1 — cxp [I3p]). 

From Eq. (29) we calculate J^^d1lna{k)^, again using Eq. (37). We find 



2 dq 
27r 



Finite temperature Drude weight of an integrable Bose chain 
and therefore 



lim Z)(r)=/3 lim r dllnA{q) 



dq 
27r 



{a+l)(3 



4cosh2 i2= 



a /j9Alna(g)f (39) 



Eq. (39) gives the high-temperature Drude weight as a function of /3/u, via 
Eqs. (38,37). Actually, we are interested in D{T) at constant particle density n, 

df_,lnaiq) dq 



-d.g = T I - 



cxp [— hia{q)] lix 

From the integral equation for 5^1na(fc) it follows, using Eq. (37), that 
dp. lna(g)|f = /3(1 + a), so that for T ^ oo, 

1 



n = o 



e-/''^ - 1 

l+n 



u(T^ oo) = - Tin 

n 

With these findings, we rewrite Eq. (39) as 



t1-^(^)-2T^ r[5.1na(.)]^^. (40) 
T-+00 2 + n + n ^ J_- 2w 



Especially, for ?^ ^ oo and rj 0, 

(1 " 
2(2+ 



lim D(T) = { 3 ' ■ . (41) 



4-2. Numerical results 

The Drude weight D{T) is calculated numerically according to Eq. (27), by setting 
up separate integral equations for the partial derivatives. The convolution integrals 
are treated using the Fast Fourier Transform algorithm. The chemical potential is 
adjusted for all temperatures considered in order to keep n at a fixed value. Results 
for g = and q = exp[— 0.5] are shown in Fig. 1. Furthermore, in Fig. 2 the next- 
leading contribution D{T) — VcK^ at low temperatures is plotted, yielding evidence 
that D{T) -VcK^ = (T^), in agreement with Eq. (24). 



5. Conclusion and outlook 



We have calculated exactly the finite-temperature Drude weight D{T) of a one- 
dimensional bosonic lattice system. Due to the conserved particle current, D{T) is 
finite. Note that the converse is not true: In the critical spin-1/2 XX Z chain, the 
current is not conserved, but D{T) is finite nevertheless. 

At low temperatures, agreement was found between the results derived from an 
effective Gaussian field theory and the BA. At high-temperatures, D{T) approaches 
a finite value which was calculated analytically from the BA solution as well. A 
numerical solution of the TBA-equations extrapolates between the low- and high- 
temperature limits. 

Open questions for future work are threefold: First, it would be interesting to 
show the absence of "Umklapp"-like terms in the g-Bose model explicitly. Next, the 
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Figure 1. The Drude weight evaluated numerically from Eq. (27), for a) g = 
and h) q = exp[— 0.5] . The curves show the results for n = 0.5, 1, 2, 4 (from bottom 
to top). The dashed lines near T = and T = 1 are the low- and high-temperature 
limits, Eqs. (22,35) and (39,41), respectively. 



amplitudes of the next-leading terms of D(T) in the low-temperature regime are of 
intrinsic theoretical interest. They are directly linked to the amplitudes of irrelevant 
operators in the effective field theory. Finally, after the particle current, the thermal 
current in the Bose chain is an interesting object for future studies. It is known 
that in other integrable models, like, for example, in the spin-1/2 XXZ-chain, the 
thermal current is conserved [19], resulting in a finite thermal Drude weight at finite 
temperatures. In how far this also holds for the q-Bose chain is to be explored in 
future work. 
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